7  V  77/ 


T_»i  nj!  r.jr  rw*/-»nr  WV  r.’T;  rijfl'^1  ■*  K.1  tj ^"TVE  V*  *7VV 


o 


om 


AD-A192  882 


HDL-CR-  88-052-2 

March  1988  f  ’■-. 


Research  in  Electromagnetic  Shielding  Theory: 
Part  2.  EMP  Simulation  Using  Small  Loops 


by  Richard  L.  Monroe 


wv 

l-r 


S 


DTIC 

ELECTE 

APR  1  9  1988 


Prepared  by 


Sol  Telecommunications  Services,  Inc. 
2000  Corporate ■ Ridge ,  Suite  810 
McLean,  VA  22102 


V 

V 

V  -* 

v 


Under  contract 

j 

DAAL02-86-C-0052 


U.S.  Army  Laboratory  Command 

Harry  Diamond  Laboratories 

Adelphi,  MD  20783*1197 

Approved  for  public  release;  distribution  unlimited. 


■«  o 

i.  <-* 


SECURITY  C.ASS  f, CATION  of  Th;S  PAGE 


REPORT  DOCUMENTATION  PAGE 


la  REPORT  SECURITY  CLASSIFICATION 

UNCLASSIFIED 
2a.  security  classification  authority 


2b.  DEClASSIFiCATlON/DOWNGRADlNG  SCHEDULE 


Form  Approved 
OMB  No  070*  0189 
Ftp  Dole  JunJO  I9SC 


ID.  RESTRICTIVE  MARKINGS 

None 

3.  DISTRIBUTION /AVAIlABiIiTy  of  REPORT 

Unlimited 


4.  PERFORM, NG  ORGANIZATION  REPORT  NUMBER(S) 


S.  MONITORING  ORGANIZATION  REPORT  NUMBERS) 

HDL-CR-88-052-2 


6a.  NAME  OF  PERFORMING  ORGANIZATION 
Sol  Telecaimunications  Services 


Cc.  ADDRESS  (Cuy,  Stale,  and  ZIP  Code) 


6b  OFFICE  SYMBOL  7a.  NAME  OF  MONITORING  ORGANIZATION 
(If  applicable) 

USALABCOM 


2000  Corporate  Ridge,  Suite  810 
Mclean,  VA  22102-7805 


7b.  ADDRESS  (City,  State,  and  ZIP  Code) 

2800  Powder  Mill  Road 
Adelphi ,  MD . 20783-1197 


Ba  NAME  OF  FUNDING  /  SPONSORING 
ORGAN, ZATiON 

USALABCOM 


£c,  ADDRESS  (Cry.  Slate,  and  ZIP  Code) 

2800  Powder  Mill  Road 
Adelphi,  MD  20783-1197 


6d  OFFICE  SYMBOL  9  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 
(If  applicable) 

AMSLC  D AAL  02-86-C-0052 


10.  SOURCE  OF  FUNDING  NUMBERS 


PROGRAM  PROJECT 

ELEMENT  NO.  NOJL162 

6.2120  120-AH25 


WORK  UN.T 
ACCESS. ON  NO 


11  TiTtc  ( Include  Security  Claudication) 

Research  in  Electromagnetic  Shielding  Theory: 

Part  2.  EMP  Simulation  Using  Small  Loops  _ 


12.  PERSONAL  AUTmOR;S) 

Richard  L.  Monroe  _ 


13a  TYPE  Of  REPORT  13b  TIME  COVERED  14  DATE  OF  REPORT  (Year,  Month,  Day)  IS.  PAGE  COUNT 


final 


i6.  Supplementary  notation 


1  3d  TIME  COVERED  1C  DATE  of  report  (Year, 

from  May  1  8  6tq  Mar  1  88  1  98  8  March  31 


II D  XE76E3;  AMSCODE:  612120.1125 


COSATl  codes 


GROUP  SUB-GROU? 


18.  SUBJECT  TERMS  (Continue  on  reverse  if  rpcessary  and  identify  by  block  number) 

E-IP  ^electromagnetic  pulse}  /  electromagnetic  shielding," 
simulation,  shielded  enclosures/  high-voltage  pulse 
generator,  small  loop  antenna,  MIL-STD-285,  IEEE-299. 


19.  ABSTRACT  (Continue  on  reverse  if  necessary  and  identify  by  block  number)  ' . 

The  feasibility  of  an  EMP  simulator  consisting  of  a  high-voltage  pulse  generator 
driving  a  small  loop  antenna  is  investigated.  It  is  shown  that  such  a  simulator 
oould  be  used  to  make  direct  time  domain  measurements  of  the  shielding  effectiveness 
of  small  enclosures  like  the  S280C  tactical  electronics  shelter  using  a  method 
analogous  to  the  frequency  domain  technique  employed  in  the  small  loop  test  described 
in  MIL-STD-285  and  IEZE-299.  Specifications  are  given  for  a  simulator  that  would 
reproduce  the  waveforn,  structure,  and  polarization  of  an  EMP  field  over  a  4  m2  plane 
surface  with  a  peak  field  of  2.7  A/m.  The  simulator  consists  of  a  2  MV  generator 
with  a  0.755  pf  storage  capacitor  and  a  resistance  of  331  ohms  driving  a  0.3  m  square 
loop  formed  from  conductors  with  1  an  radii.  In  operation,  the  loop  is  oriented  with 
its  plane  perpendicular  to  the  surface  of  the  enclosure  and  its  center  2.5  m  from 
that  surface. 


20.  DISTRIBUTION /AVAILABILITY  OF  ABSTRACT 
0  UNCI  ASS  IF  IE  DAJ  Nil  Ml  TED  □  SAME  AS  RPT 
22a  NAME  OF  RESPONSIBLE  INDIVIDUAL 

W.  Coburn 


DD  FORM  1473,  84  MAR  83APRe 


[TcT  I  21.  ABSTRACT  SECURITY  CLASSIFICATION 

:  as  rpt.  □  d tic  users _ same  as  _ 

22b  TELEPHONE  (include  Area  Code)  22c.  OFFICE  SYMBOL 

(703)  490-2320  SLCHD-NW-EC 


83  APR  edit, on  may  be  used  until  exhausted  security  Classification  Of  this  page 

All  other  editions  are  obsolete 


1 


CONTENTS 


A 


1.  INTRODUCTION, 


2.  MAGNETIC  FIELD  DISTRIBUTION  OF  A  SMALL  LOOP, 


3.  PULSE  GENERATION  IN  AN  INDUCTIVE  LOAD, 


4.  SPECIFICATIONS  FOR  A  HEMP  SIMULATOR, 


5.  CONCLUSIONS, 


LITERATURE  CITED, 


DISTRIBUTION 


Acotsic'.  To- 


NUS  CPA&i 
DUO  TAE 
U  '  O  i‘" 

JiiSiif'Cviiioii 


i  DijiiiL-uti'  n  / 


I  Avj  hiL.lOy  CcJt-3  “1 


W'/ 


»•»* ;  or 

~  j| 


FIGURES 


m v. v r  i." *■  n*v  »v 


1.  The  S280C  tactical  electronics  shelter . . . 

2.  Average  HEMP  field:  (a)  frequency  domain  and 

(b)  time  domain . . . 

3.  A  rectangular  loop  composed  of  a  conductor  with  radius 

r  and  linear  dimensions  a  and  b  lying  in  the  y,z  plane 
of  a  rectangular /cyl indrical  coordinate  system . 

4.  Magnetic  field  versus  frequency  at  the  point 

x  =  y  =  z  =  0  due  to  a  uniform  harmonic  current  of 
unit  magnitude  on  a  loop  antenna  with  dimensions 
a  =  b  =  0.3  m  lying  in  the  y,z  plane  with  its  center 
located  a  distance  d  =  0.48  m  from  the  x,y  plane  . 

5.  Magnetic  field  distribution  in  the  x,y  plane  due  to  a 
uniform  harmonic  current  of  unit  magnitude  on  a  loop 
antenna  with  a  =  b  =  0.3  m  lying  in  the  y,z  plane  with 
its  center  located  a  distance  d  =  1  m  from  the  x,y  plane. 


page 

8 


i 


12 


4 


15  4 


16 


18  ! 

i 

t 


v  •. 


\y.  ■- 


£  , 

a  -*»  -S  .V  -VA-»  -V 


6.  Magnetic  field  distribution  in  the  x,y  plane  due  to  a 


uniform  harmonic  current  of  unit  magnitude  on  a  loop 


antenna  with  a  =  b  =  0.3  m  lying  in  the  y,z  plane  with 


its  center  located  a  distance  d  =  1.5m  from  the  x,y  plane 


7.  Magnetic  field  distribution  in  the  x,y  plane  due  to  a 


uniform  harmonic  current  of  unit  magnitude  on  a  loop 


antenna  with  a  =  b  =  0.3  m  lying  in  the  y,z  plane  with 


its  center  located  a  distance  d  =  2  m  from  the  x,y  plane. 


8.  Magnetic  field  distribution  in  the  x,y  plane  due  to  a 


uniform  harmonic  current  of  unit  magnitude  on  a  loop 


antenna  with  a  =  b  =  0.3m  lying  in  the  y,z  plane  with 


its  center  located  a  distance  d  =  2.5  m  from  the  x,y  plane 


9.  Magnetic  field  distribution  in  the  x,y  plane  due  to  a 


uniform  harmonic  current  of  unit  magnitude  on  a  loop 


antenna  with  a  =  b  =  0.3  m  lying  in  the  y,z  plane  with 


its  center  located  a  distance  d  =  3  m  from  the  x,y  plane 


10.  Schematic  (a)  and  equivalent  circuit  (b)  of  a  pulse 


generator . 


11.  Load  current  I  produced  by  a  2  MV  pulse  generator 


where  Cg  =  0.755  pf ,  R  =  331  n,  and  L  =  0.69  ^h 


WW 


12.  Load  currents  I  produced  by  a  2  MV  pulse  generator 


where  L  =  0.69  /ih  ,  Cg  =  0.1  nf . 

13.  Load  currents  I  produced  by  a  2  MV  pulse  generator 

where  L  =  0.69  /ih ,  Cg  =  1  nf . 

14.  Load  currents  I  produced  by  a  2  MV  pulse  generator 
=  0.69  /Jh,  C£  = 


where  L 


10  nf 


INTRODUCTION 


1  . 

A  recent  study*  has  suggested  that  small  loops  driven  by  pulse 

generators  could  be  used  to  create  transient  source  fields  for  the 

purpose  of  measuring  the  efficiency  of  small  tactical  shelters  like 

the  S280C  (fig.  1)  as  shields  against  EMP  (electromagnetic  pulse) 

fields.  This  source,  together  with  appropriate  sensors  at  designated 

locations  inside  the  enclosure,  would  comprise  a  time  domain  analog 

2  3 

of  the  small  loop  test  described  in  MIL-STD-285  and  IEEE-299.  As 
prescribed  by  these  standards,  the  loop  would  be  located  close  to 
the  shelter  with  its  plane  perpendicular  to  the  plane  of  the  shelter 
wall  in  order  to  maximize  coupling  to  the  interior.  The  objective  of 
the  pulsed  loop  would  be  to  generate  a  field  on  a  designated  area  of 
the  enclosure  surface  that  matches  the  structure  and  waveshape  of 
the  EMP  so  that  the  efficiency  of  that  region  as  an  EMP  shield  can 

*R.  L.  Monroe,  A  Theory  of  Electromagnetic  Shielding  with 
Applications  to  MIL-STD  285,  IEEE-299,  and  EMP  Simulation,  Harry 
Diamond  Laboratories,  HDL-CR-85-052- 1 ,  (February  1985). 

2 

Department  of  Defense,  MIL-STD-285,  Method  of  Attenuation 
Measurements  for  Enclosures,  Electromagnetic  Shielding,  for 
Electronic  Test  Purposes  (25  June  1956). 

3 

IEEE,  Proposed  IEEE  Recommended  Practices  for  Measurement  of 
Shielding  Effectiveness  of  High  Performance  Shielding  Enclosures, 
IEEE  299,  IEEE  Inc.,  NY  (June  1969). 


be  determined  directly.  With  this  approach,  the  overall  efficiency 
of  the  enclosure  would  be  determined  by  repeating  the  measurements 
at  different  locations  until  all  surfaces  exposed  to  the  BMP  are 
covered . 

The  use  of  a  small  loop  as  a  piece-wise  EMP  simulator  in  this 
manner  is  possible  because  the  electric  and  magnetic  fields  inside  a 
shielded  enclosure  produced  by  any  external  source  depend  only  on 
the  tangential  components  of  the  source  magnetic  field  on  the  out¬ 
side  surface  of  the  enclosure^.  That  is,  the  internal  fields  are 
virtually  independent  of  the  source  electric  field.  Consequently, 
any  two  sources  that  produce  the  same  or  nearly  the  same  tangential 
magnetic  fields  on  the  outside  of  the  enclosure  will  be  indistin¬ 
guishable  from  inside  the  enclosure,  and  any  source  that  reproduces 
another  source's  tangential  magnetic  field  on  the  surface  of  the 
enclosure  can  simulate  that  source  as  far  as  internal  electric  and 
magnetic  fields  are  concerned.  Thus  an  efficient  magnetic  field 
source  like  the  loop  can  simulate  an  EMP  over  a  designated  region  by 
approximating  the  structure  and  waveshape  of  the  EMP ’ s  tangential 
magnetic  field  without  any  need  to  simulate  its  electric  field.  A 
comparison  of  the  internal  fields  with  the  simulated  EMP  magnetic 
field  from  this  source  can  then  be  taken  as  a  measure  of  the 
effectiveness  of  that  region  as  a  shield  against  EMP  fields. 

*R.  L.  Monroe,  A  Theory  of  Electromagnetic  Shielding  with 
Applications  to  MIL-STD  285,  IEEE-299,  and  EMP  Simulation,  Harry 
Diamond  Laboratories,  HDL-CR-85-052-1 ,  (February  1985). 


The  decoupling  of  internal  electric  and  magnetic  fields  from 
external  electric  fields  occurs  because  shielded  enclosures  are 
composed  of  materials  belonging  to  the  class  of  good  conductors  and 
are  constructed  in  such  a  way  that  all  openings  (nonconducting  areas 
in  the  enclosure  walls)  are  incorporated  in  structures  that  are 
equivalent  to  waveguides  beyond  cutoff.  These  measures  insure  that 
the  net  electric  field  (incident  plus  reflected  field)  at  the  sur¬ 
face  of  the  enclosure  is  very  small  and  is  determined  not  by  the 
electric  field  of  the  external  source  but  by  an  impedance  function 
that  characterizes  the  surface  of  the  enclosure.  The  relationship 
between  the  impedance  function  and  the  electric  and  magnetic  fields 
at  the  surface  is  expressed  by  impedance  boundary  conditions. 
Shielded  enclosures  must  satisfy  such  conditions  at  all  points  in 
order  to  function  properly.  If  an  enclosure  should  fail  to  satisfy 
impedance  boundary  conditions  because  of  damage  or  poor  design,  this 
failure  will  be  revealed  by  severely  degraded  shielding  performance. 

In  this  note,  we  will  investigate  the  feasibility  of  small 
loops  as  simulators  of  high-altitude  EMP  (HEMP)  fields  over  a  por¬ 
tion  of  an  enclosure  surface.  The  characteristics  of  the  HEMP 
magnetic  field  that  the  simulator  must  reproduce  are  its  polariza¬ 
tion  (direction),  spatial  variation  (change  in  magnitude  as  a  func¬ 
tion  of  position),  and  waveshape  (time  history).  Since  HEMP  fields 
incident  at  the  surface  of  the  earth  can  be  approximated  as  uniform 
plane  waves  with  both  vertical  and  horizontal  components  over  scale 


10 


sizes  corresponding  to  the  dimensions  of  the  S280C  and  with  average 

4 

time  histories  like  that  shown  in  figure  2  ,  the  basic  function  of 
the  simulator  must  be  to  generate  a  spatially  uniform,  unipolar, 

transient  magnetic  field  in  any  desired  direction  with  a  risetime 

_  9 

(time  to  peak  field)  of  10  ns  =  10  x  10  s  or  less.  For  purposes  of 
testing  the  shielding  properties  of  enclosures,  it  is  not  necessary 
that  the  simulator  produce  peak  magnetic  fields  as  large  as  the 
HEMP,  although  such  a  capability  would  be  very  useful  for  other 
types  of  tests.  The  basic  objective  here  is  to  produce  a  large 
enough  external  field  with  the  correct  waveform  over  a  sufficiently 
large  area  of  the  enclosure  so  that  the  internal  fields  can  be 
measured . 

In  the  following  section,  we  will  show  that  a  small  loop  can 

indeed  produce  a  uniform  magnetic  field  over  a  significant  area  on  a 

nearby  plane  surface  and  that  the  direction  of  the  field  can  easily 

be  varied  from  horizontal  to  vertical.  This  demonstration  will 

reveal  the  basic  tradeoff  between  field  strength  and  exposure  area 

that  controls  the  performance  of  this  type  of  simulator.  Section  3 

will  examine  circuit  paarmeters  for  a  pulse  generator  to  drive  the 

loop  with  a  transient  current  that  will  produce  a  magnetic  field 

with  the  EMP  waveshape.  The  results  of  these  sections  will  then  be 

combined  in  section  4  to  give  specifications  for  a  device  that  will 

2 

simulate  the  HEMP  magnetic  field  over  an  area  of  4  m  with  a  peak 
field  of  at  least  1  A/m. 


4 

Bell  Laboratories,  EMP  Engineering  and  Design  Principles,  Whippany 
NJ  (1975). 
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Figure  2.  Average  EMP  field:  (a)  frequency  domain  and  (b)  time 


domain  (reproduced  from  ref.  4). 


MAGNETIC  FIELD  DISTRIBUTION  OF  A  SMALL  LOOP 


2  . 


A  loop  is  considered  small  as  a  radiating  structure  if  its 
enclosed  area  A  is  such  that 

A1/2  <  <  A  =  C  /f  (2.1) 

o 

where  A  is  the  wavelength  of  its  transmitted  field,  C  is  the  velo- 

o 

city  of  light,  and  f  is  the  frequency  of  the  source.  When  this  con¬ 
dition  is  satisfied,  the  current  distribution  on  the  loop  is  uniform 
( independennt  of  position)  and  its  fields  differ  little  from  those 
of  a  point  dipole  source  even  at  locations  in  the  extreme  near 
field. *  Since  A  is  inversely  proportional  to  f,  (2.1)  can  be  used  to 
estimate  the  maximum  allowable  size  for  any  small  loop  over  a  spec¬ 
ified  range  of  frequencies.  Thus  a  small  loop  simulator  for  the  HEMP 
field  shown  in  figure  2  will  be  limited  by 

A1/2  <  <  3  x  108/1  xlO8  =  3  m  (2.2) 

since  this  field  includes  a  continuous  range  of  frequencies  from  0 

g 

to  10  Hz.  The  preceding  gives 

a  a  0 . 3  m  (2.3) 

as  an  upper  limit  for  the  size  of  a  small  square  loop  over  the 
frequency  range  of  interest.  In  the  following,  we  will  observe  the 
limit  imposed  by  (2.3)  and  consider  only  small  square  loops  less 


*R.  L.  Monroe,  A  Theory  of  Electromagnetic  Shielding  with 
Applications  to  MIL-STD  285,  IEEE-299,  and  EMP  Simulation,  Harry 
Diamond  Laboratories,  HDL-CR-85-052-1 ,  (February  1985). 

13 


than  or  equal  to  0.3  m  on  a  side  as  candidates  for  the  simulator. 
One  could,  of  course,  consider  larger  loops  or  loops  with  different 
shapes.  However,  it  is  believed  that  these  would  offer  no  advant¬ 
ages.  In  fact,  larger  loops  would  inevitably  increase  the  inductance 
in  the  driver  circuit-  a  result  that  is  usually  undesirable.  Thus 
the  natural  tendency  is  to  consider  smaller  rather  than  larger 
loops.  As  far  as  shape  is  concerned,  it  appears  to  have  little  ef¬ 
fect  on  the  performance. 

For  design  purposes,  the  small  square  loop  offers  the  additional 
advantage  that  its  near-field  distribution  is  relatively  easy  to 
compute  using  expressions  developed  in  Monroe^  (equ.  (5.39)).  These 
expressions  are  referred  to  the  rectangular/cylindrical  coordinate 
shown  in  figure  3  where  the  plane  of  the  loop  lies  in  the  y,z  plane 
with  its  center  C(0,0,d)  on  the  z  axis  a  distance  d  from  the  origin. 
With  this  arrangement,  the  field  distribution  in  the  x,y  plane  due 
to  a  uniform  harmonic  current  with  amplitude  I  on  the  loop  can  be 
computed  by  setting  z  =  0  in  (5.39)  and  then  specifying  the  frequen¬ 
cy  f,  the  dimensions  of  the  loop  a  =  b,  its  distance  from  the  x,y 
plane  d  =  H,  and  the  coordinates  of  the  point  in  the  x,y  plane  where 
the  field  is  to  be  computed.  Figure  4  is  a  plot  of  the  magnitude 
H(0,0)  of  the  magnetic  field  at  the  origin  versus  f  computed  from 


*R.  L.  Monroe,  A  Theory  of  Electromagnetic  Shielding  with 
Applications  to  MIL-STD  285,  IEEE-299,  and  EMP  Simulation,  Harry 
Diamond  Laboratories,  HDL-CR-85-052-1 ,  (February  1985). 
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Figure  3.  A  rectangular  loop  composed  of  a  conductor  with  radius  r 


and  linear  dimensions  a  and  b  lying  in  the  y,z  plane  of  a 
rectangular/cyl indr ical  coordinate  system. 


I  "V  _  »»  iJTf 


Figure  4,  Magnetic  field  versus  frequency  at  the  point  x=y=z=0  due  to  a 
uniform  harmonic  current  of  unit  magnitude  on  a  loop  antenna  with 
dimensions  a  =  b  =  .3m  lying  in  the  y,z  plane  with  its  center  located 
a  distance  d  =  .48  m  from  the  x,y  plane. 


16 


4 


equation  (5.39)  for  a  harmonic  current  of  unit  amplitude  (I 


A)where  a  =  b  =  0.3  m  and  d  =  0.48  m.  This  plot  illustrates  another 


useful  property  of  the  small  loop,  namely,  a  flat  response  over  a 


very  wide  range  of  frequencies  which  in  this  case  includes  the  EMP 


frequencies  0  to  1  x  10  Hz.  This  means  that  when  the  harmonic  cur¬ 


rent  source  is  replaced  on  the  loop  by  a  transient  source  with  a 


bandwidth  of  100  MHz  or  less  the  waveshape  of  the  magnetic  field  in 


the  x,y  plane  will  be  an  image  of  the  current  on  the  loop.  With  this 


loop,  one  can  adjust  the  parameters  of  a  driver  circuit  to  produce  a 


transient  current  with  the  correct  EMP  waveform  and  be  assured  that 


the  current  will  generate  a  magnetic  field  in  the  x,y  plane  with  the 


same  waveform.  We  will  take  advantage  of  this  fact  in  the  following 


section  to  specify  a  pulse  generator  for  the  simulator. 


The  flat  response  of  the  loop  also  means  that  we  can  compute 


the  magnetic  field  distribution  throughout  the  entire  x,y  plane 


using  (5.39)  at  a  single  frequency,  say  1  MHz,  and  that  distribution 


will  be  valid  for  all  frequencies  in  the  range  0  to  100  MHz.  This 


has  been  done  to  obtain  figures  5  through  9,  where  the  continuous 


lines  forming  three  nested  rectangles  are  scale  drawings  of  the 


side,  front,  and  door  of  the  S280C  shelter.  The  line  segments  origi¬ 


nating  on  a  rectangular  grid  of  points  spaced  0.2  m  apart  represent, 


by  their  length  and  direction,  the  magnitude  and  direction  of  the 


magnetic  field  in  the  x,y  plane  at  those  points.  In  this  representa¬ 


tion,  the  magnitude  of  the  field  at  each  point  is  normalized  to  the 


field  H ( 0 , 0 )  at  the  origin  (0,0)  which  is  assigned  a  length  equal  to 


the  grid  spacing  0.2  m.  Since  H(0,0)  is  always  the  maximum  field  in 


"«*  j,’*  ,*•  M**  y*  J 
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w. 


luS 
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f=  1  MHz 


d3  1m 


Tick  marks  are  spaced  .2  n  apart  on 
both  axes 


H(0,G)a  7.40  x  10  A/m 


Projection  oP  loop  on  x,y  plane 


6.9  x  10 


Figure  5.  Magnetic  field  distribution  in  the  x,y  plane  due  to  a  uniform 


harmonic  current  of  unit  magnitude  on  a  loop  antenna  with  a  s  b  s  .3m 


lying  in  the  y,z  plane  with  its  center  located  a  distance  d  s  1  m  from 


the  x,y  plane. 


Figure  6.  Magnetic  field  distribution  in  the  x,y  plane  due  to  a  uniform 
harmonic  current  of  unit  magnitude  on  a  loop  antenna  with  a  s  b  s  ,3  o 
lying  in  the  y,z  plane  with  its  center  located  a  distance  d  s  1,5  m  from 
the  x,y  plane. 
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Figure  7.  Magnetic  field  distribution  in  the  x,y  plane  due  to  a  uniform 


harmonic  current  of  unit  magnitude  on  a  loop  antenna  with  a  s  b  =  .3m 


lying  in  the  y,z  plane  with  its  center  located  a  distance  d  =  2  m  from 


the  x,y  plane. 


V  %  .•//*  lV  .» 


aw 


f=  1  MHz 


d=  3m 

HCBjQ)®  2,66  x  10  4  A/m 


6.9  X  10 


Tick  narks  are  spaced  .2  n  apart  on 
Jgj/DDth  axes 


Projection  oP  loop  on  x,y  plane 


Figure  9.  Magnetic  field  distribution  in  the  x,y  plane  due  to  a  uniform 
harmonic  current  of  unit  magnitude  on  a  loop  antenna  with  a  =  b  =  .3m 
lying  in  the  y,z  plane  with  its  center  located  a  distance  d  =  3  m  from 
the  x,y  plane. 
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the  x,y  plane  for  this  arrangement,  the  fields  at  other  points  are 
always  represented  by  lines  that  are  less  than  or  equal  to  0.2  m  in 
length.  Thus,  a  field  in  the  x  direction  that  is  uniform  in  magni- 
itude  over  a  certain  area  will  be  seen  here  as  a  series  of  lines 
spaced  .2  m  apart  parallel  to  the  x  axis  running  through  that  area. 

Figure  5  shows  the  field  distribution  in  the  x,y  plane  from  a 
square  loop  0.3  m  on  a  side  located  in  the  y,z  plane  with  its  center 

on  the  z  axis  1  m  from  the  origin.  In  addition  to  the  field  distri¬ 

bution,  the  figure  displays  I,  f,  d,  H(0,0),  a,  and  b,  as  defined 
previously.  Also  shown  are  r,  the  radius  of  the  loop  conductor,  and 
L,  the  loop  inductance,  where  the  latter  is  computed  using  the  form- 

5 

ula  found  in  Grover.  In  this  case,  a  1  A  current  on  the  loop  gives 

_3 

H(0,0)  =  7.4  x  10  A/m  which  means  that  a  pulse  generator  producing 
a  peak  current  of  1  kA  =  1000  A  on  the  loop  would  produce  a  magnetic 
field  in  the  x  direction  with  a  peak  of  7.4  A/m  at  (0,0).  Unfortun¬ 
ately,  the  figure  also  indicates  that  the  area  over  which  the  field 

2 

can  be  described  as  uniform  is  less  than  1  m  .  The  latter  is  consid¬ 
erably  smaller  than  an  earlier  estimate*  for  this  setup  which  was 


5 

F.  W.  Grover,  Inductance  Calculations:  Working  Formulas  and 
Tables,  Dover  Publications,  NY  (1962). 

*  R.  L.  Monroe,  A  Theory  of  Electromagnetic  Shielding  with 
Applications  to  MIL-STD  285,  IEEE-299,  and  EMP  Simulation,  Harry 
Diamond  Laboratories,  HDL-CR-85-052- 1 ,  (February  1985). 
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based  on  calculations  at  an  insufficent  number  of  locations. 


To  improve  the  field  di3tribution ,  it  is  necessary  to  move  the 


loop  farther  from  the  x,y  plane.  This  has  been  done  in  figure  6 


where  d  has  been  set  equal  to  1.5  m  without  changing  the  other  input 


parameters.  The  figure  does  show  a  significant  improvement  with  a 


uniform  field  that  covers  virtually  the  entire  area  of  the  door  (1.3 


m  ).  However,  this  improvement  has  been  obtained  at  the  cost  of  a 


reduction  in  H(0,0)  from  7.4  x  10  ^  to  2.2  x  10  ^  A/m.  Thus,  for  a 


given  loop  current,  there  is  a  tradeoff  between  the  area  that  can  be 


covered  with  a  uniform  field  and  the  magnitude  of  that  field.  Fig¬ 


ures  7,  8,  and  9  show  how  still  larger  areas  of  uniform  field  can  be 


obtained  by  moving  the  loop  farther  from  the  x,y  plane.  Of  these, 


figure  8  where  d  =  2.5  m  appears  to  satisfy  our  objective  of  a  uni¬ 


form  field  over  an  area  of  4  in  .  Since  a  1  A  loop  current  gives 


H ( 0 , 0  )  =  4.6  x  10  A/ra  in  this  case,  a  peak  1/op  current  of  2.2  kA 


will  be  required  to  produce  a  simulated  EMP  with  the  desired  peak 


field  of  1  A/m , 


These  figures  indicate  that  the  field  in  the  x,y  plane  has  some 


of  the  characteristics  of  a  magnetic  dipole  with  virtual  poles  lo¬ 


cated  on  the  x  axis  equidistant  from  the  y  axis.  The  area  of  uniform 


field  lies  between  the  poles  which  move  farther  apart  as  the  loop 


moves  away  from  the  the  x,y  plane.  The  field  between  the  poles  is 


predominately  in  the  x  direction  Thus,  if  the  positive  y  axis  is 


defined  as  pointing  in  the  vertical  direction,  the  setup  described 


here  corresponds  to  a  vertically  polarized  EMP  field  (vertical 


■  '  V  V  V  V.  V  V 


electric  field,  horizontal  magnetic  field)  incident  on  the  side  or 


front  of  the  shelter.  Other  polarizations  can  be  obtained  simply  by 


rotating  the  plane  of  the  loop  about  the  z  axis  by  the  appropriate 


number  of  degrees.  A  90  rotation  puts  the  plane  of  the  loop  in  the 


x,z  plane  and  produces  a  uniform  magnetic  field  in  the  y  direction 


corresponding  to  a  horizontally  polarized  EMP  field. 


v;/,v  v, 


PULSE  GENERATION  IN  AN  INDUCTIVE  LOAD 


The  preceding  has  shown  that  it  will  be  necessary  to  generate  a 


transient  loop  current  with  a  peak  of  at  least  several  thousand  amps 


on  a  time  scale  of  10  ns  in  order  to  satisfy  the  simulator  require¬ 


ments  we  have  specified.  This  must  be  done  by  driving  the  loop  with 


a  high-voltage  pulse  generator  such  as  the  one  described  schematic¬ 


ally  in  figure  10  (a).  The  basic  operation  of  these  devices  is  well 


known  and  will  not  be  discussed  here  beyond  noting  that  the  process 


is  initiated  by  charging  capacitor  Cg  to  a  voltage  V  using  the  rec¬ 


tifier  circuit  to  the  left  of  the  switch  S.  When  Cg  is  fully  charg¬ 


ed,  S  is  thrown  to  its  lower  position,  discharging  Cg  to  ground  and 


generating  a  current  pulse  in  the  load  impedance  Z.  An  equivalent 


circuit  giving  an  approximate  (linear)  description  of  the  current 


pulse  is  shown  in  figure  10  (b)  where  V  is  considered  to  be  applied 


instantaneously  at  t  =  0  on  the  input  terminals.  In  this  circuit,  R. 


and  L.  refer  to  the  internal  resistance  and  inductance  of  the 
l 


pulser,  and  R  and  L  refer  to  the  external  resistance  and  induct- 
r  e  e 


ance  of  the  load.  These  can  be  lumped  together  by  defining  the  total 


resistance  and  inductance  of  the  circuit  as  follows: 


R .  =  R .  +  R 

tie 


(3.1) 


L.  =  L.  +  L 
tie 


For  a  purely  inductive  load  like  the  small  loop,  we  can  expect  R^  >> 


R  and  L  >>  L.  so  that  (3.1)  can  be  rewritten 
e  e  i 


»n.rtsPw,r»»fc 


where  R  is  resistance  of  the  pulser  and  L  is  inductance  of  the  loop, 


Solving  the  circuit  for  the  current  gives 


I ( jw  )  = 


(3.3) 


jw  [R  +  juL  +  1/jwC  ] 


6 

The  latter  is  easily  transformed  to  the  time  domain,  giving: 


For  D  >  0 


I(t)  =  I i exp ( -C j t/2 ) s inh ( D  t) 


(3.4) 


For  D  <  0 


I ( t )  =  I^expi -C^t/2 ) sin[ ( -D )  t] 


(3.5) 


where 


Ix  =  V/[LD1/2) 

I2  =  V/[L(-D)1/2] 

D  =  (C./2)2  -  C 
1  o 

C  =  1/(LC  ) 
o  s 

Cx  =  R/L  . 


(3.6) 

(3.7) 

(3.8) 

(3.9) 


(3.10) 


Thus  I ( t )  will  have  either  an  exponential  or  a  damped  sinusoidal 
waveform  depending  on  whether  Cg,  R,  and  L  are  such  that  D  is  great¬ 
er  than  or  less  than  zero.  Since  the  HEMP  waveform  is  itself  of  the 


G.  Doetsch,  Guide  to  the  Application  of  Laplace  Transforms,  D.  Van 
Nostrand  Company  LTD,  London,  (1961). 


exponential  type,  our  interest  here  will  be  limited  to  the  former 
case,  and  we  will  attempt  to  choose  the  circuit  parameters  so  that 
I(t)  approximates  the  curve  in  figure  2.  In  doing  this,  we  will  have 
great  freedom  in  choosing  the  internal  parameters  Cg  and  R.  Both  of 
these  parameters  can  be  varied  over  several  orders  of  magnitude  if 
necessary.  However,  our  choice  of  L  is  much  more  limited.  We  can 
always  increase  L  by  adding  a  series  inductor  or  using  multiturn 
loops.  However,  this  is  not  a  useful  option,  since  it  would  increase 
the  circuit  impedance  and  decrease  the  current.  We  can  decrease  L  by 
decreasing  the  size  of  the  loop  while  increasing  the  diameter  of  the 
loop  conductor.  However,  size  reductions  in  high-voltage  systems  are 
limited  by  electrical  breakdown,  and  it  appears  unlikely  that  an  in¬ 
ductance  much  below  the  value  for  the  0.3  m  square  loop,  L  =  6.9  x 

-7 

10  h  =  0.69  ph  could  be  realized  in  a  practical  simulator.  Accord¬ 
ingly,  we  will  set  the  simulator  specifications  using  L  =  0.69  ph  as 
the  target  inductance  for  the  system. 


SPECIFICATIONS  FOR  A  HEMP  SIMULATOR 


Since  the  HEMP  waveform  in  figure  2  can  be  expressed  analytic- 

__  _ i _ a. _  _  _  •  .  _ _ _ _ *  ■  i  4 


ally  as  the  difference  between  two  decaying  exponentials, 
H  =  H  [°xp(-a.t)  -  exp(-a  tl]  A/m 

O  1  L 

where 


(4.1) 


H  =  5.25  x  10/377  =  139  A/m 
o 


a1  =  4  x  106  s"1 


4.76  x  108  s-1 , 


(4.2) 

(4.3) 

(4.4) 


it  is  useful  to  write  equation  (3.4)  in  the  same  form: 


I(t)  =  I ^exp ( -C i t/2 ) sinh { D  t 


=  I  exp(-C  t/2 ) [exp(D1/2t)  -  exp ( -D1 /2t ) ] / 2 


=  (Ij/2)  exp[-(C1/2-D1/2)tJ  -  exp[-(C1/2+D1/2)t]  .  (4.5) 


We  can  then  compare  (4.1)  and  (4.5)  and  note  that  these  expressions 
will  describe  the  same  waveform  provided 


al  " 

Cj/2 

- 

D1/2 

a2  = 

Cl/2 

+ 

D1/2 

where  D 

and 

C 

^  are 

are 

f ixed , 

(4. 

6) 

can 

ing 

C  and 

.  R. 

If 

(4.6 

s 

will 

speci 

f  y  a 

pulser 

(4.6) 


Bell  Laboratories,  EMP  Engineering  and  Design  Principles,  Whippany 
NJ  (1975). 
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t' 


that  reproduces  the  HEMP  waveshape.  Since  (4.6)  is  nonlinear  m  R 


and  Cs>  we  must  resort  to  successive  approximations  to  obtain  a 


solution.  Using  an  initial  estimate  of  R  =  300  ft  and  C  =  .5  nf ,  we 

s 


obtain 


R  =  331  ft  , 


C  =  0.755  nf  . 


(4.7) 


The  reader  can  verify  that  (4.7)  satisfies  (4.6)  to  an  accuracy  of 


better  than  1  %, 


Figure  11  is  a  plot  of  I(t)  computed  with  (3.4)  for  a  pulse 


generator  specified  by 


V  :  2  X  106  V  :  2  MV 


R  =  331  ft 


(4.8) 


C  =  .755  nf 
— s - 


driving  an  inductive  load  equal  to  that  of  a  0.3  m  square  loop 


where 


L  =  0.69  fj h 


(4.9) 


It  shows  a  peak  current  of  5.8  kA  and  a  waveform  identical  to  the 


HEMP  field  in  figure  2.  This  pulser  driving  the  0.3  m  square  loop 


would  generate  a  spatially  uniform  magnetic  field  with  the  HEMP 


waveform  and  a  peak  field  of  2.7  A/m  over  an  area  of  4  m  on  a 


plane  surface  located  2.5  m  from  the  center  of  the  loop  (fig.  8). 


Thus,  it  would  satisfy  the  requirements  for  the  EMP  simulator  speci¬ 


fied  at  the  outset. 


„»  v  7S  •  s  'j  v  v  v  v  'j  '.-  v  ■j.'.-V  ■  j>„.  v  ".-A-  v'/v 


t  (nanoseconds) 


Figure  11  .  Load  current  I  produced  l>y  a  2  MV  pulse  generator 
where  C  =  .755  pf,  R  =  331  n,  and  L  =  .69  u h. 


CONCLUSIONS 


5  . 


The  simulator  we  have  specified  would  allow  one  to  determine 

the  overall  efficiency  of  an  enclosure  as  a  shield  against  HEMP 

2 

fields  by  a  piece-wise  process  in  which  4  m  sections  of  the  enclo¬ 
sure  are  exposed  in  succession  to  magnetic  fields  with  a  HEMP-like 
structure,  polarization,  and  waveform.  The  peak  external  field  of 

2.7  A/m  should  produce  a  dynamic  range  of  at  least  65  dB  for  shield- 

-  3 

ing  measurements  since  an  internal  field  of  1  x  10  A/m  =  1  mA/m 
can  usually  be  measured.  If  greater  dynamic  range  is  needed,  one  can 
obtain  this  at  the  cost  of  a  reduced  area  of  coverage  by  moving  the 
loop  closer  to  the  enclosure.  For  example,  by  moving  the  center  of 
the  loop  from  d  =  2.5  m  to  d  =  1.5  m  (fig.  6),  the  peak  field  is  in¬ 
creased  by  a  factor  of  4.7  to  12.7  A/m.  which  translates  to  a  dynam¬ 
ic  range  of  more  than  80  dB .  The  shielding  effectiveness  of  the  vast 

majority  of  shelters  will  be  found  to  fall  within  this  range. 

Finally,  we  note  that  it  may  be  useful  to  simulate  HEMP  fields 
with  waveforms  differing  from  the  waveform  shown  in  figure  2  since 
the  latter  is  only  an  average.  This  can  be  done  easily  by  adjusting 

C  and  R.  However  there  are  limitations  on  what  can  be  achieved  in 
s 

this  way.  Figures  12,  13,  and  14  show  the  types  of  waveforms  that 

can  be  obtained  by  letting  C  range  from  0.1  to  10  nf  and  R  from  200 

s 

to  1600  ft.  In  general,  these  figures  show  that  smaller  R  and  larger 

Cg  produces  larger  peak  currents  occurring  at  later  times  while 

larger  R  and  smaller  C  produces  smaller  peaks  at  earlier  times 


compared  to  figure  11.  Thus,  an  earlier  (<  10  ns)  peak  can  be 
obtained  by  sacrificing  magnitude,  and  a  larger  magnitude  can  be 
obtained  by  accepting  a  later  {>  10  ns)  peak,  but  a  peak  that  is 
both  larger  and  earlier  than  the  one  in  figure  11  is  not  possible 
unless  L  can  be  significantly  reduced. 


t  (nanoseconds) 


Figure  12.  Load  i-urrents  1  produced  bv  a  2  M\  pulse  generator 
where  I.  r  .09  ph  ,  C  ~  .1  nf,  and: 

* 1 .  R  =  200  n 

* 2 .  R  =  400  fi 

*3 .  R  =  800  n 

#4  R  =  1000  n 
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t  (nanoseconds) 

F igure 

13.  Load  currents 

I  produced  by  a  2  MY  pulse  generator 

where 

L  =  .  6  9  u  h ,  f  =] 

S 

n  f ,  and 

#1  . 

R  =  200  n 

#2 

R  =  300  a 

*3 

R  =  400  n 

*4 

R  =  800  a 

I  (Rfl) 


1  10  100 
t  (nanoseconds) 

Figure  14.  Load  currents  7  produced  by  a  2  MV  pulse  generator 
where  L  =  .69  yh ,  -  10  nf,  and 

# i  R  =  200  n 

#2  R  =  300  a 

#3  R  =  4  00  ft 


#4  R  =  800  ft 
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ATTN:  TECHNICAL  DIRECTOR,  AMSLC-TD 

2800  Powder  Mill  Road 
Adelphi,  MD  20783-1145 


Installation  Support  Activity 

ATTN:  Library,  SLCIS-IM-TL  (3  copies) 

ATTN:  Library,  SLCIS-IM-TL  (Woodbridge) 

ATTN:  Legal  Office,  SLCIS-CC 

2800  Powder  Mill  Road 

Adelphi,  MD  20783-1145 

USAISC 


ATTN:  Record  Copy,  ASNC-ADL-TS 


ATTN:  Technical  Reports  Branch,  SLCIS-IM-TR 

2800  Powder  Mill  Road 
Adelphi,  MD  20783-1145 


45 


A  v  v%  A  A  ■  %.  %  , 


Harry  Diamond  Laboratories 


ATTN:  D/Division  Directors 


2800  Powder  Mill  Road 
Adel phi ,  MD  20783-1197 

Harry  Diamond  Laboratories 
ATTN:  Chief,  SLCHD-NW-E 

ATTN:  Chief,  SLCHD-NW-EC 

ATTN:  Chief,  SLCHD-NW-ED 

ATTN:  Chief,  SLCHD-NW-EE 

ATTN:  SLCHD-NW-EC,  W.  Coburn,  COTR  (9  copies) 

ATTN:  SLCHD-NW-RE ,  A.  Cuneo 

2800  Powder  Mill  Road 
Adelphi ,  MD  20783-1197 


